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ABSTRACT 

In this letter we present ground-based subarcsecond mid- infrared imaging and spectroscopy of young 
super star clusters in the overlap region of the merging galaxies NGC4038/4039 (the Antennae) 
obtained with the VLT Imager and Spectrometer for mid-Infrared (VISIR). With its unprecedented 
spatial resolution VISIR begins to resolve the Hn/PDR complexes around the star-forming regions 
for the first time. In the N-band spectra of two young star clusters unexpectedly low polycyclic 
aromatic hydrocarbon (PAH) emission is observed, compared to what is seen with the Infrared Space 
Observatory (ISO) and with the Spitzer Space Telescope. We conclude that a large fraction of the PAH 
emission cannot directly be associated with the super star clusters, but originate from an extended 
region of at least 215 pc radius around the clusters. In the distribution of [Neil] 12.81 fim emission a 
highly obscured cluster is discovered that does not have an optical or near- infrared counterpart. 
Subject headings: galaxies: individual (NGC4038/4039) — galaxies: star clusters — ISM: dust/PAH 
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1. INTRODUCTION 

Starburst galaxies experience a phase of rapid evolu- 
tion. The rate at which their gas reservoir is turned 
into stars cannot be maintained for long, which makes 
the starburst phase by definition a transient one. The 
resulting stellar clusters make starbursts unique labora- 
tories for the study of star formation, stellar populations 
and the evolution of galaxies as a whole. 

Since the earliest stages of massive star formation are 
generally heavily enshrouded by dust, infrared observa- 
tions are essential to detect the youngest stellar popu- 
lations and provide various diagnostic features to study 
the properties of the interstellar matter (ISM) and the 
underlying stellar population. Recently, a new genera- 
tion of ground- and space-based instruments working at 
mid-infrared wavelengths has become operational, giving 
this field of research a large impulse. 

The Antennae (NGC4038/4039, A rp244) is the near- 
est m ajor merger of two large spirals ijToomre fc Toomrd 
1972). Since the beginning of the interaction the sys- 
tem went through several episodes of violent star for- 
mation of which the last one is probably still ongoing 
ijVigroux et alJll99fil) . The resulting (super) star clusters 
have been studied e xtensively throughout the electro- 
magn etic spectrum llWhitmore et all 120051 iWang et all 
120041 iGilbert et al] 12000ft . Radio and mid-infrared ob- 
servations showed that the region between the two spi- 
rals (the overlap region) hosts spectacular obscured star 
formation. The brightest mid-infrared component pro- 
duces 15% of the total 15 fim luminosity of the entire 
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syste m HMhahel et al .1119981 IHummel fc Van der Hulstl 
1986). This region is covered by a prominent dustlanc, 
and can be associate d with a faint, red source in th e HST 
images (number 80 of Whitmore fc Schweizerll995l here- 
after WS95), illustrating how optical data alone are in- 
sufficient to identify and study the youngest star-forming 
regions. An ISOCAM spectrum of this region shows 
strong fine-structure emission lines ([Neil] and [Nein]) 
and pronounced em ission from polycycl ic aromatic hy- 
drocarbons (PAHs) <|Mirabel et al]ll998jh 

We use the recently commissioned VL T Imager and 
Spect rometer for mid-Infrared (VISIR) l)Lagage et alJ 
120041) at the Very Large Telescope (VLT) of the Euro- 
pean Southern Observatory (ESO) to study the most lu- 
minous super star clusters and the surrounding matter 
in the Antennae overlap region in detail. 

2. OBSERVATIONS AND DATA REDUCTION 

The southern part of the overlap region of 
NGC4038/39 was observed with VISIR on 2005 
April 17 - 20. VISIR offers a large set of imaging 
filters covering the 8-13 and the 17 - 24 micron 
range of which we used four (see Table all at the 
largest pixel scale of O'.'127/pixel, resulting in a 32.5"x 
32.5" field-of-view. To correct for the strong and highly 
variable atmospheric background chopping and nodding 
was applied, with a chopper throw of 14". The total 
on-source integration time was 60 minutes in both 
the PAH and PAH reference filter, 120 minutes in the 
[Neil] filter and 15 minutes in the [Neil] reference filter. 
The latter observations had to be aborted prematurely 
because of strong winds. Between individual nod cycles 
the telescope was given a random jitter offset to correct 
for array artefacts. Stand ard stars (HD99167 and 
HD157999. fCohen et al.lll999ft were observed before and 
after each target observation for flux calibration. 

Additionally, low-resolution N-band spectra were ob- 
tained for the two brightest sources in the overlap region 
simultaneously (R « 185 - 390, 0"75 slit). Four spectral 
settings overlapping by at least 15% were used to cover 
the full N-band (A contr = 8.8, 9.8, 11.4 and 12.4 fim). 
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Fig. 1 — Left panel: SoFI Ks-band (2.16 (im). The VISIR slit is 
shown in white, north is up, east is left. Right panel: VISIR [Nell] 
12.81 fim of the overlap region, inserted: source la and lb in the 
[Neil] 12.81 jim filter with the contours of PAH 11.25 /im filter 
emission overlayed. 



TABLE 1 

Cluster photometry and sizes 



filter 


-^ccntr 


AA 


seeing 8 


flux density 13 


size c 




(jjm) 


(/an) 




(mjy) 


(PC) 








Source 1 






PAH 


11.25 


0.59 


0'.'33 ± tf'03 


98 ± 20 


88 ± 14 


PAHref 


11.88 


0.37 


0'.'32 ± tf'01 


133 ± 27 


80 ± 14 


[Nenjref 1 


12.27 


0.21 


0'.'46 ± O'.'Ol 


203 ± 41 




[Neil] 


12.81 


0.18 


0^34 ± 0'.'02 


373 ± 75 c 


93 ± 20 








Source 2 






PAH 


11.25 


0.59 


0'/33 ± 0^'03 


74 ± 15 


76 ± 9 


PAHref 


11.88 


0.37 


0'.'32 ± O'.'Ol 


84 ± 17 


71 ± 9 


[Nenjref 1 


12.27 


0.21 


0'.'46 ± O'.'Ol 


114 ± 23 




[Neil] 


12.81 


0.18 


0'/34 ± 0^'02 


210 ±42 


84 ± 10 



a mean of FWHM of PSF of all standard stars observed before 
and after target observations 

b total flux density, sum of compact and broad component 

C FWHM of the broad component of the best fit 

d low SNR; observations aborted prematurely because of strong 
wind 

c sum of sources la and lb, one-fifth can be associated with lb 



The total on-source integration time was 30 minutes for 
each spectroscopic setting. Since the images showed no 
large-scale diffuse emission, chopping and nodding on-slit 
was applied with a 10" chopper throw. Early-type stars 
were observed before and after each target observation 
for flux calibration (HD115892 and HD135742). 

IRAF 1 routines and customized IDL scripts were used 
for data reduction. Subtracting the chopping and nod- 
ding pairs removes most of the sky and telescope back- 
ground. The resulting images were shifted and coad- 
ded. Wavelength calibration and removal of the curva- 
ture in the spectroscopic data was done by tracing sky- 
lines. Spectra were extracted in a l'/27 aperture. Ab- 
solute flux calibration was obtained by normalizing the 
standard star spectra to VISIR narrow-band fluxes. The 
accuracy of the flux calibration is 20% both for imaging 

1 Image Reduction and Analysis Facility (NOAO, National Op- 
tical Astronomy Observatories) 
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Fig. 2. — Normalized radial profiles of source 1 (upper panels) 
and 2 (lower panels) in the PAH (left panels) and reference contin- 
uum filters (right panels). Filled circles show the profile of source 
1/2, open circles that of all standard stars combined. Error bars 
are plotted, but are smaller than the plotting symbol, expect for 
the points at smallest radial distance. 
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Fig. 3. — Upper panel: low resolution N-band spectrum of source 
la. The apparently enhanced noise from 9-10 fim is a result of 
the log scale of this plot. The observed fine-structure line ratio 
[Slv]/[Nell] is higher than obtained within the larger apertures 
of ISO and Spitzer, indicating that VISIR, probes regions of higher 
excitation closing in on the location of most intense star formation. 
Lower panel: idem source 2. 
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and spectroscopy. 

3. RESULTS 

The narrow-band PAH 11.25 fim and reference im- 
ages show two bright resolved sources separated by 5'.' 8 
(w 600 pc in projection at a distance of 21 Mpc 2 , as- 
suming a Hubble constant of 70 km s" 1 Mpc -1 ; at 
this distance 1" corresponds to 105 pc) in the south- 
ern part of the overlap region. These sources correspond 
to the two stro ng peaks identified in the overlap region 
in ISO images l)Mirabel et al]ll998[K The eastern source 
is the counterpart of the faint, very red cluster 80 in 
HST images (WS95) and a bright cluste r in the near- 
infrared (IDwmc 157. lBrandl etafl l2005). The western 
source can be associated with a very crowded, much less 
obscured, bright cluster complex consisting of the star 
clusters WS95-86 to WS95-90 plus a number of close 
neighbours and with another bright near-infrared clus- 
ter, IDwirc 136. It is remarkable how similar the east- 
ern and western source are in shape and brightness in the 
mid-infrared, while these regions appear very different at 
shorther wavelengths. Given the similar ages for these 
cluste rs (both < 6 Mvr lGilbert et al.ll2000HMengel et all 
2005), the differences in appearance at the near-infrared 
and the optical indicate strongly varying physical condi- 
tions and/or extinction in these star-forming regions. 

In the [Ne n] filter image the eastern source breaks up 
into two distinct sources, which are separated by Of! 54 
(« 60 pc in projection, Fig. ^ source la & lb). Rela- 
tive astrometry strongly suggests that source la can be 
associated with the bright near-infrared cluster IDwirc 
157 and with the optical source WS95-80 mentioned be- 
fore. This was determined by matching source 2 through 
all wavelength bands: the peak in the mid-infrared data 
was accurately matched with the well-defined peak in the 
near-infrared. The morphology of the cluster complex in 
the near-infrared was in turn precisely overlayed with 
the optical data. The fainter second source lb, emitting 
one-fifth of the total emission in the [Ne n] filter (contin- 
uum plus [Neil] line flux), has not been observed before. 
Neither in the optical nor in the near-infrared possible 
counterparts are identified. At shorter wavelengths this 
cluster is hidden by extinction and at longer wavelengths 
most instruments lack spatial resolution to separate the 
two clusters. 

Under the assumption that source lb is as la younger 
than 6 Myr, that source lb has a similar spectral en- 
ergy distribution and using the values for the Ks-band 
magnitude and extinction as published (Ks = 14.59 and 
Ay is 4.23, Mengel et al. 2005; assuming that Ay/Ax 
w 9.3. IMathislll99Cij) . the minimal extinction necessary 
to account for the observations is A K > 7.8 (Ay > 72) 
of source lb. Although the exact intrinsic Ks-band flux 
of source lb is not straightforward to determine, it is 
obvious that it is extinguished by an enormous amount 
of obscuring matter. The presence of this cluster com- 
promises efforts to model the emission from this region 
with one coeval stellar population. It might for example 
partly explain the discrepancy between the mid-infrared 
luminosity measured in this region by ISO (15% of the 

2 Note that the distance of the Antennae is und er debate. Re- 
cently a lower distance of 13.8 ± 1.7 Mpc was found ISaviane et alJ 
2004), which would affect the values derived here. 



TABLE 2 
Line fluxes from spectroscopy 



— ; 

Species 




'Vest 




line flux 


— — 

EW 




(fim) 


/-in — 14 — 1 — 9 \ 

(10 erg s 1 cm 


(cm) 






Source la 




[Ar III] 


8.99138 


3.9 ± 0.9 


0.033 ± 0.009 


[Siv] 


10.5105 


7.6 ± 0.9 


0.057 ± 0.011 


PAH 


11.3 


14.0 ± 3.8 


0.057 ± 0.017 


[Neil] 


12.8136 


49.6 ± 6.2 


0.092 ± 0.018 






Source 2 




[Ar m] 


8.99138 


4.2 ± 0.7 


0.051 ± 0.011 


[Siv] 


10.5105 


10.3 ± 2.3 


0.078 ± 0.019 


PAH 


11.3 


25.0 ± 4.2 


0.127 ± 0.029 


[Neil] 


12.8136 


30.2 ± 4.5 


0.111 ± 0.023 



total 15^m flux of the entire system, iMirabel et al.ll99^) 
and the far lower expected mid-infrared luminosity as re- 
constructed from the properties of the underlying stellar 
population derived from nea r-infrared and optic al data 
(3.6% of the total 15/im flux. IMengellt~aTll2005]) . 

Fig. |21 shows that both source 1 and 2 are spatially re- 
solved. The radial profiles of the sources are well fitted 
by a combination of two Gaussians, one compact compo- 
nent and a more extended component, both contributing 
roughly 50% to the total flux. The full- width at half max- 
imum (FWHM) of the broad component for the best fit 
varies from 0'.'85 - l'.'O for source 1 and slightly smaller, 
0'.'75 - 0'.'90, for source 2. This corresponds to a size 
of 90 - 105 pc and 80 - 95 pc respectively, similar to 
sizes observed for giant molecular clouds. Table ^ lists 
the size of the broad component and flux density F„ for 
both sources in all filters. 

The spatial distribution along the slit in the spectro- 
scopic data shows two bright, slightly resolved sources 
without any obvious more extended emission. The ex- 
tracted N-band spectra show the typical signatures of 
a region of massive star formation: fine structure lines 
([Nell], [Siv] and faint [Arm]) and a rising contin- 
uum characteristic for thermal emission by dust. Near- 
infrared data suggest that these clusters are very young, 
which is supported by the presence of strong [S iv] emis- 
sion. Line fluxes and equivalent widths (EWs) are given 
in Tabled 

Besides obvious similarities in the general shape of the 
N-band spectra (Fig.^J) there are some remarkable differ- 
ences from pr evious observations o f the same regions by 
ISO (Fig.l in IMirabel et alJll998j) and Spitzer (Brandl 
et al, in preparation). The most striking contrast is 
the strength of the PAH features. Where the ISO and 
Spitzer spectra of the region around complex 1 show a 
pronounced 11.3 fim PAH- feature, the VISIR spectrum 
only shows a slight indication of it. 

4. PAH EMISSION 

One of us (BRB) is leading a GTO program to ob- 
serve the Antennae with Spitzer-IRS. The results will be 
published in a subsequent paper (Brandl et al. in prepa- 
ration) . These observations include observations of both 
complex 1 and source 2. 

Comparison of the VISIR spectrum of source la pub- 
lished here with the Spitzer spectrum of the same region 



4 



Snijders et al. 



shows considerable differences. The continuum flux den- 
sity measured by VISIR is roughly three-quarters of that 
observed by Spitzer. Of the total emission in the 11.3 
fim PAH feature VISIR measures a much smaller frac- 
tion, only approximately one-quarter of the emission seen 
in the Spitzer spectrum. This large difference indicates 
that the region of PAH emission is significantly more ex- 
tended than that of continuum emission. 

Because of its wider slit Spitzer is expected to collect 
more flux than VISIR in case of extended sources (4'.' 7 
and 0'.'75 respectively). Indeed, VISIR finds a lower flux 
density than Spitzer at all wavelengths, but the relative 
fraction is different for the various spectral features. The 
EW of the 11.3 ^m PAH feature found with VISIR is over 
three times smaller than the EW measu red within the 
large r apertures of both Spitzer and ISO l)Mirabel et alJ 
1998). The total flux in the PAH feature in the VISIR 
spectrum of source la is (1.4 ± 0.4) ■ 10~ 13 ergs -1 cm~ 2 . 
The remaining three-quarters of PAH emission that is 
"missing" compared to the Spitzer data originate from 
a more extended region. If we assume that this (4.2 
± 0.8) • 10 -13 erg s _1 cm~ 2 of PAH emission is ho- 
mogeneously distributed and has a surface brightness 
just below the detection limit of our data (2a is 2.56 
• 10 -14 erg s _1 cm~ 2 arcsec~ 2 ), the emission must orig- 
inate from a region of at least 215 ± 25 pc radius. 

While this faint emission cannot directly be detected 
in our data, indications for more extended PAH emission 
are also found by close examination of the radial profiles 
in both the PAH and the reference continuum filter where 
a slightly softer drop-off of the profile in the PAH filter 
is found (see Fig. 0J. 

These results indicate that an extended PAH emission 
component below the VISIR surface brightness detection 
limit is present that accounts for the discrepancy of the 
VISIR with the ISO and Spitzer data. Spitzer IRAC im- 
ages indeed show in dications of widespread PAH emission 
l|Wang et al.ll2004() . 

For source 2 there are no ISO data available, but com- 
paring the VISIR and Spitzer data leads in a similar way 
to the conclusion that extended PAH emission must be 
present. The VISIR spectrum of this source shows a more 
pronounced PAH feature. The EW is roughly two times 
smaller than that measured in the Spitzer spectrum and 
of the total PAH flux one-third is measured by VISIR. 
Compared to source la we measure a smaller fraction of 
the total continuum flux, only two-thirds of the contin- 
uum flux density measured by Spitzer. The smaller frac- 
tion of the continuum measured and the smaller discrep- 
ancy between the VISIR and Spitzer EWs reflects the 
extended structure of the underlying rich cluster com- 
plex. In contrast to the case with the relatively simple 
morphology of source 1, here a smaller fraction of the spa- 
tially extended structure is covered by the narrow VISIR 
slit. Although the morphology of source 2 is simple in 
our imaging data, the fact that we only measure half of 
the continuum compared to Spitzer indicates that a more 
extended component must be present below the VISIR 
detection limit. Taking this into account the PAH emis- 
sion must still originate from a more extended region 
than the continuum. 

These results indicate that VISIR is resolving the 
Hn/PDR complexes. While the ionized gas and con- 
tinuum emission, and thus the hard radiation field of the 
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Fig. 4. — Normalized radial profile of source 1 in the PAH (filled 
circles) and reference continuum niters (open circles). The pro- 
file in the PAH filter drops off slightly softer than the continuum 
profile. Error bars are plotted, but are smaller than the plotting 
symbol, expect for the points at smallest radial distance. 



most recent star formation, are much more concentrated, 
a large fraction of the PAH emission originates from a 
region more extended than 215 pc around the clusters. 
This could be explained by the presence of a (slightly) 
more evolved, less massive and more widespread pop- 
ulation of field stars by which the PAH molecules are 
excited. This population can consist of the original, pre- 
merger disk stars or can be the result of dissolution of 
clusters formed in previous episodes of star formation. 
The timescale of the latter process is estimated to be a 
few 10 Myr in the Antennae (|Mengel et al JI2005D . In this 
scenario a large fraction of the PAH emission is not di- 
rectly associated with the currently observed super star 
clusters and therefore does not trace the very most re- 
cent star formation. PAH emission has been detected i n 
less UV-rich enviro nments before ljUchida et alJll998j) . 
iLi fc Draind l)2002j) argue that it does not necessarily 
take hard UV radiation to excite PAHs and produce the 
observed emission bands in the 6-12 /jm range. Another 
possible explanation is the destruction of PAH molecules 
by hard UV photons up to large distances from the ioniz- 
ing clusters, w hich is also claimed to be seen in other star- 
burst galaxies l|Tacconi-Garman et alJ200 5: Beira c~e"t alJ 
2006). However, source 2 shows the strongest PAH emis- 
sion and has the highest [Siv]/[Neii] ratio (0.15 ± 0.02 
and 0.34 ± 0.04 for sources la and 2 respectively). This 
significantly higher ratio, indicative of a harder radiation 
field, would argue against the PAH destruction scenario 
in this case. 

5. FUTURE WORK 

The VISIR N-band observations of the Antennae over- 
lap region presented in this letter show the power of 
a ground-based mid-infrared instrument mounted on a 
large telescope. Due to the unprecedented spatial res- 
olution in the mid-infrared, which space observatories 
lack, VISIR is an excellent tool to use in detailed studies 
of nearby starburst galaxies. Structures on the scale of 
giant molecular clouds (80 - 105 pc) are resolved, each 
hosting one or a number of star clusters. 

In order to further investigate the excitation of PAHs 
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in starburst environments, we are currently studying a 
number of starburst galaxies with a range in age and 
luminosity. 



We thank the Paranal Observatory Team for their sup- 
port. 



REFERENCES 



Beirao, P., Brandl, B. R., Dcvost, D., Smith, J. D., Hao, L. & 

Houck, J. R., 2006, ApJ 643, 1 
Brandl, B. R. et al. 2005, ApJ 635, 280 

Cohen, M., Walker, R. C, Carter, B., Hammcrsley, P., Kidger, M. 

& Noguchi, K., 1999, AJ 117, 1864 
Gilbert, A. M. et al. 2000, ApJ 533, L57 
Hummel, E. & Van der Hulst, J. M., 1986, A&A 155, 151 
Lagage, P. O. et al. 2004, ESO Messenger No. 117, 12 
Li, A. & Draine, B. T., 2002, ApJ 572, 232 
Mathis, J. S., 1990, ARA&A 28, 37 

Mengel, S., Lehnert, M. D., Thatte, N. & Genzel, R., 2005, A&A 
443, 41 

Mirabel, I. F. et al. 1998, A&A 333, LI 

Neff, S. G. & Ulvestad, J. S., 2000, AJ 120, 670 



Rubin, V. C, Ford, W. K. J. & D'Odorico, S., 1970, ApJ 160, 801 
Saviane, I., Hibbard, J. E. & Rich, M. R., 2004, AJ 127, 660 
Tacconi-Garman, L. E., Sturm, E., Lehnert, M., Lutz, D., Davics, 

R. I & Moorwood, A. F. M., 2005, A&A 432, 91 
Toomre, A. & Toomre, J., 1972, ApJ 178, 623 
Uchida, K. I., Sellgren, K. & Werner, M., 1998, ApJ 493, 109 
Vigroux, L. et al. 1996, A&A 315, L93 
Wang, Z. et al. 2004, ApJ 154, 193 
Whitmore, B. C. & Schweizer, F., 1995, AJ 109, 960 
Whitmore, B. C, Zhang, Q., Leitherer, C, Fall, M. S., Schweizer, 

F. & Miller, B. W., 1999, AJ 118, 1551 
Whitmore, B. C. et al. 2005, AJ 130, 2104 



